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ABSTRACT 

In this research, a methodology to predict damage initiation, damage growth, fatigue 
life and residual strength in titanium matrix composites (TMC) is outlined. Emphasis was 
placed on micromechanics-based engineering approaches. Damage initiation was predicted 
using a local effective strain approach. A finite element analysis verified the prevailing 
assumptions made in the formulation of this model. Damage growth, namely, fiber-bridged 
matrix crack growth, was evaluated using a fiber bridging (FB) model which accounts for 
thermal residual stresses. This model combines continuum fracture mechanics and 
micromechanics analyses yielding stress-intensity factor solutions for fiber-bridged matrix 
cracks. It is assumed in the FB model that fibers in the wake of the matrix crack are idealized 
as a closure pressure, and an unknown constant frictional shear stress is assumed to act along 
the debond length of the bridging fibers. This frictional shear stress was used as a curve 
fitting parameter to the available experimental data. Fatigue life and post-fatigue residual 
strength were predicted based on the axial stress in the first intact 0° fiber calculated using the 
FB model and a three-dimensional finite element analysis. 
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NOMENCLATURE 


Current crack length and unbridged initial crack length, m 

Hole diameter, m 

Crack growth coefficient constants 

Fiber and matrix modulus, MPa 

Composite longitudinal, transverse and shear modulus, MPa 
Boundary correction factor for applied stress intensity factor 

Weight function 

Strain concentration factor 

Orthotropic stress concentration factor 

Discrete stress intensity factor range in matrix, MPa yjm 

Continuum stress intensity factor range in composite, MPa's/m 


Distance between parallel cracks, m 
Range in bridging pressure, MPa 
Fiber radius, m 
Hole radius, m 

# omin / 

Composite stress ratio = 0 / Qmn 


Matrix stress ratio = 



Applied stress range, MPa 

Change in temperature, °C 
Fiber volume fractions 


Width of specimen, m 

Coefficient of thermal expansion for fiber and matrix, mm/mm/ °C 
Effective strain parameter 

Strain range 

Interfacial frictional shear stress, MPa 

Composite Poisson’s ratio 

Axial stress in fiber and matrix, MPa 

Axial stress in the fiber at and , respectively, MPa 

Failure stress of first intact fiber under static loading, MPa 

Failure stress range in first intact fiber under fatigue loading, MPa 

Axial thermal residual stresses in the matrix, MPa 
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INTRODUCTION 


Titanium matrix composites (TMC) have excellent potential as structural materials for 
elevated temperature applications such as future hypersonic flight vehicles. However, a major 
problem facing TMC is the lack of a generalized life prediction methodology that can be 
applied in designs incorporating these advanced materials. The complex damage growth 
process that evolves in these materials can limit the use of conventional damage tolerance 
design approaches often used for metals. On the composite microstructural level, damage 
consists of multiple failure mechanisms including fiber-matrix debonding, fiber breakage and 
matrix cracking. The interaction and combination of these failure mechanisms can lead to 
extensive damage process zones having no single dominate crack. In such instances, classical 
fracture mechanics approaches are not suitable. An additional complexity are thermal residual 
stresses that are induced in the fiber and matrix during the fabrication process which can have 
a profound effect on the fracture and fatigue behaviors of TMC. The heterogeneous nature of 
TMC poses challenges to the development of a damage tolerance design philosophy for this 
class of materials. 

In an effort to address these challenges, a micromechanics-based life prediction 
methodology is proposed. The effects of applied stress, number of fatigue cycles and 
specimen geometry (e.g. stress concentrations and notches) on damage initiation, damage 
growth and residual strength can be determined and could be used as a design limiter for 
notched TMC. Figure 1 is a flow chart describing this methodology and how it can be used to 
predict the following fatigue scenarios: (1) initiation of cracks from notches using a local 
effective strain parameter; (2) growth of matrix cracks to threshold condition where crack 
arrest occurs based on stress intensity factor in the matrix calculated using a fiber bridging 
(FB) model; (3) continuous growth of matrix cracks based on stress intensity factor in the 
matrix using the FB model; (4) specimen life based on the stress in the first intact fiber 
(unmachined fiber nearest the notch or hole) calculated using the FB model; and (5) post- 
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fatigue residual static strength (static strength after fatigue loading) based on the fiber stress 
calculated using either the FB model or a three-dimensional finite element analysis. 

The objectives of this research is to outline the framework for this generalized life 
prediction methodology. This paper will review research addressing issues of damage 
initiation, damage growth and post-fatigue residual strength. Emphasis was placed on 
micromechanics-based approaches. Damage initiation was predicted using a local effective 
strain approach developed by Hillberry and Johnson [1]. A finite element analysis verified the 
prevailing assumptions made in the formulation of this model. Next, damage growth and life 
prediction approaches for fiber bridged matrix crack growth were developed based on a fiber 
bridging (FB) model. The FB model combines continuum fracture mechanics and 
micromechanics analyses to derive stress intensity factor solutions of matrix cracks bridged by 
fibers. Thermal residual stresses were accounted for in this model by Bakuckas and Johnson 
[2]. The frictional shear stress that is assumed to act along the debond length of the bridging 
fibers was used as a curve fitting parameter for the available experimental data. Finally, post- 
fatigue residual static strength was predicted based on the axial stress in the first intact 0° fiber 
calculated using the FB model developed in [2] and using a three dimensional finite element 
analysis of Bigelow et al. [3,4,5] 

EXPERIMENTAL PROCEDURES 
Materials and Specimens 

Two TMC systems were studied in this research, namely, SCS-6/TM5-3 and SCS- 
6/Timetal-21S laminates. Both material systems were fabricated by Textron by hot isostatic 
pressing (HIPing) titanium foils between unidirectional tapes of silicon-carbide (SCS-6) fibers 
having a diameter of 0.14 mm. In the SCS-6/Ti-15-3 and SCS-6/Timetal-21S laminates, 
fibers were held in place with molybdenum and titanium-niobium cross weave wires, 
respectively. The composition of matrix materials used in the SCS-6/Ti-15-3 and SCS- 
6/Timetal-21S laminates, were Ti-15V-3Cr-3Al-3Sn and Ti-15Mo-3Al-2.7Nb-0.2Si, 
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respectively. For the SCS-6/TM5-3 laminates, the lay-ups studied were [0]g and [0/90] s with 
a range of measured fiber volume fractions from 0.33 to 0.36. For the SCS-6/Timetal-21S 
laminates, the lay-ups studied were [0/90] s and [0/90]2 S with a range of measured fiber volume 
fractions v f from 0.36 to 0.38. The material properties are listed in Table 1. 

Specimens were machined using a diamond wheel saw into straight-sided coupons with 
the 0° fibers in the loading direction. Each specimen was 152.5-mm in length. Center holes 
having diameters (d ) of 6.3±0.3-mm were drilled ultrasonically in specimens having a width of 
19.1-mm. Center slits having lengths (2a) of 5.08-mm, 7.62-mm and 8.89-mm and a width of 
0.254-mm were machined via electo-discharged machining (EDM) in specimens having a 
width of 25.4-mm. The SCS-6/Ti-15-3 specimens were tested in the as-received condition 
while the SCS-6/Timetal-21S specimens were heat treated at a temperature of 620°C for 8 
hours in a vacuum and then fumaced cooled prior to testing. To enhance optical observations 
and surface replication, the surface of each specimen was polished to obtain a flat and lustrous 
finish. Thin sheet aluminum end tabs having dimensions 19.1-mm x 19.1-mm were bonded 
on the ends of all specimens in order to prevent grip failure. Information on the specimens 
tested is listed in Table 2. 


Test Procedures 

Testing was conducted using closed-loop servo-hydraulic test machine equipped with 

water cooled hydraulic grips. Tension-tension fatigue tests were done under load control with 
R Comp = 0.1 ( 5nu % miX ) at a frequency of 10 Hz at both room temperature and an elevated 

temperature of 200°C. Results from these tests were generated in previous studies [3,6,7]. In 
these studies, damage initiation and progression were monitored and recorded using several 
techniques including a closed-circuit television system (CCTV), a long focal length 
microscopic system, surface and edge replicas, scanning electron microscopy (SEM), and 
optical microscopy. In addition, the acoustic emission (AE) technique was used where 
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emphasis was placed on detecting the occurrence of fiber breakage. Details of the AE set-up 
and parameters are reported in [8]. 

ANALYTICAL PROCEDURES 

Figure 1 is a flow chart outlining the micromechanics-based life prediction 
methodology developed. It was used in this study to determine the effects of applied stress, 
number of fatigue cycles and specimen geometry (e.g. stress concentrations and notches) on 
damage initiation, damage growth, fatigue life and residual strength. The method consists of 
three analysis approaches: (1) local strain model used to predict damage initiation based on 
the value of the effective strain parameter; (2) fiber bridging model used to evaluated damage 
growth based on the stress intensity factor in the matrix and used to evaluate fatigue life and 
post-fatigue residual strength based on the stress in the first intact fiber, and; (3) a three- 
dimensional finite element model used to predict the virgin static and post-fatigue residual 
static strengths based on the stress in the first intact fiber. These three approaches were 
developed in previous studies and were combined in this paper to make up the framework for a 
life prediction methodology. A brief description of these three approaches and the application 
to damage initiation, damage growth, fatigue life, and residual strength are discussed next. 


Matrix Crack Initiation 

A procedure to predict the initiation of damage in the form the matrix cracks was 
introduced in [1]. In the approach, an effective strain parameter was used which includes 
thermal residual stresses and stress concentrations: 


= 


K.e^ + 


c fl 


"htrm \ 


K,Ae 


( 1 ) 


■'m ) 


where K ( , £ max , , # m , and Ae are the orthotropic stress concentration factor, the 

maximum applied strain, the axial matrix thermal residual stress in the 0° plies, the matrix 
modulus, and the applied strain range, respectively. It is assumed in this model that the strain 
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in the fiber, matrix, and composite are compatible and that the strain concentration factor in 
the matrix, K e t , is equal to the orthotropic (composite) stress concentration factor, K ( . 

The axial matrix thermal residual stress, (?* m , induced during the fabrication process 
was calculated using the VISCOPLY code [9], a micromechanics based model. The thermo- 
viscoplastic properties of the constituents used in the calculation of are reported in [9] 
for SCS-6/Ti-15-3 laminates and in [10] for SCS-6/Timetal-21S laminates. In order to include 
thermal residual stresses in the calculations, the composite fabrication process was analytically 
simulated by modeling the cool down cycle from an assumed stress-free temperature to room 
temperature. A stress-free temperature of 555 °C was used for the SCS-6/Ti-15-3 laminates 
and 620°C was used for the SCS-6/Timetal-21S laminates. For SCS-6/Ti-15-3 laminates, it 
was assumed that the stress-free condition occurs at one half the melting temperature of the 
titanium matrix, that is, at 555 °C. For the SCS-6/Timetal-21S laminates, it was assumed that 
the residual stresses were removed during the heat treatment cycle at a temperature of 620°C. 
Values of a^ trm induced during cool-down of the composite fabrication are listed in Table 1. 

The effective strain parameter as a function of cycles to failure (final fracture) for the 
Ti-15-3 matrix material was determined in [1] and is shown in Figure 2. The solid line is a 
curve fit to the experimental data. For a notched composite with matrix thermal residual 
stresses, the number of cycles to matrix crack initiation can be predicted by comparing the 
effective strain parameter defined by Equation (1) with the solid line in Figure 2. 

Finite Element Analysis 

A finite element analysis was conducted in this research to verify the prevailing 
assumptions made in the effective strain approach. The assumptions made in this model 
(Equation 1) are the following: (1) strain in the fiber, matrix and composite are compatible, 
(2) the orthotropic stress concentration factor is the same as the strain concentration factor for 
the fiber, matrix and composite. Thus, the stress concentration calculated using a anisotropic 
continuum analysis is assumed to be equivalent to the strain concentration in the matrix. A 
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concern with this assumption is that the local stress concentrations that occurs near the ends of 
discrete machined fibers, particularly for hole configurations, are not accounted for. In some 
instances, cracks have been observed to initiate and grow as shown in Figure 3 where the 
initiation site is in the region of a machined fiber end. In order to determine the severity of 
the stress concentration factor around discrete fiber ends along the hole boundary, a two- 
dimensional finite element analysis was conducted for a unidirectional SCS-6/Timetal-21S 
monolayer having a fiber volume fraction, \j = 0.36, and containing a center hole, d/W = 
0.33. 

The finite element model used is shown in Figure 4 which consists of a global region 
and a local heterogeneous region. The global region was used to model the specimen 
geometry to which the far-field loads were applied and was treated as an orthotropic elastic 
continuum whose properties were calculated from the constituent properties using the rule-of- 
mixtures. The local heterogeneous region was used to conduct the micromechanics analysis 
near the center hole. In this region, the fiber and matrix was modeled as distinct zones using 
their respective neat properties. The constitutive properties used in the finite element analysis 
of the SCS-6/Timetal-21S monolayer are given in Table 1. 

The entire finite mesh consisted of 2588 elements and 10,581 nodes. Nine-noded 
isoparametric elements with 3x3 Gauss numerical integration was used. Thermal residual 
stresses were included in the analysis by simulating the cool-down cycle from an assumed 
stress-free temperature of 620°C to room temperature, as described previously. Then a 
mechanical load was applied incrementally up to 700 MPa along the top row of nodes. 

The distribution of the Von-Mises effective stresses in the matrix is shown in the stress 
contours in Figure 5. As shown in this figure, stress concentrations do exist in the matrix in 
the vicinity of the machined fiber ends. For the sake of clarity, not all contours are labeled for 
the matrix material between the first and second cut fibers and the second and third cut fibers. 
A detailed stress contour for the matrix between the first (back) cut fiber and second (front) cut 
fiber is shown in Figure 6 indicates that the stress builds-up behind the second (front) cut 
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fiber. Similar results were obtained for the matrix material between the remaining cut fibers, 
Figure 5; the stress builds-up behind the front cut fibers. Despite these local stress 
concentrations, the maximum stress and strain concentration still occurs at Location A in 
Figure 5. Thus, the center hole is the primary stress concentration while the machined fiber 
ends are secondary stress concentrations. 

The local strain concentration factor in the matrix is defined as: 


p local- max 
__ C axial 
* ~ global 
^ axial 


( 2 ) 


where is the maximum local axial strain (at Location A, Figure 5) and e s J£f is the 

remote global axial strain. For an elastic analysis from the initial load increment, Kf = 
3.81. The global stress concentration factor for a hole in an orthotropic continuum is defined 
[11] as: 



The terms in this equation are the laminate properties as defined in the nomenclature section. 
The term a is the finite width correction factor. In this study, the finite width correction 
factor for isotropic material was used [12]: 



For the composite analyzed, K, — 3.83 which is in excellent agreement with the strain 

concentration factor (K,) calculated from the discrete finite element analysis and, therefore, 
can be used in the definition of the local effective strain given by Equation (1). 
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Matrix Crack Growth 


Methods to evaluate matrix crack growth was introduced in several studies using fiber 
bridging (FB) models [13-15]. The application of these FB models to analyze damage growth 
in unidirectional TMC was demonstrated in [6]. The frictional (fiber-matrix interface) shear 
stress term assumed in these models was used as a fitting parameter to crack growth data. It 
was shown in this study that the FB models do provide an engineering tool for interpreting 
crack growth data in this class of materials. However, the values of the frictional shear stress 
term used to evaluate crack growth data varied with several parameters including the applied 
load and crack length. Further developments of fiber bridging analysis was realized when 
thermal residual stresses were considered [2]. It was shown in [2] that a single value of 
frictional shear stress was sufficient to evaluate the crack growth data for a given material 
system, lay-up and applied temperature. While these studies are encouraging, the FB models 
still lack true predictive capabilities because the frictional shear stress term must be determined 
experimentally. 

Fiber bridged matrix crack growth was analyzed using the approach outlined in [2]. In 
that study, crack growth behavior was evaluated according to a version of the power law 
function introduced by Walker which includes R-ratio effects [16]: 

f v 


da 

dN 


= C- 


AK. 


(1 -K.) 


(l - m) 


( 5 ) 


where C, n, and m are material constants for the neat matrix material, R m is the local stress 

_irun / 

ratio in the matrix (R m = 'A,, ) and A K m is the Mode I stress intensity factor range in the 


matrix. By using Equation (5), thermal residual stresses can be included in the analysis 
through the quantities R m and AK m as described next. 
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Calculation of R m 

The stress ratio in the matrix (R m = 

code [9]. After the composite fabrication simulation, a global mechanical load from S'™ to 
S™** was applied, and the corresponding stresses in the matrix in the 0° plies were calculated. 
Using these values of stress, the value of R„ in the matrix in the 0° plies was calculated. 
Values of R m are listed in Table 2 for each specimen tested. As revealed by the results in 
Table 2, the values of R„ were considerably higher than the value of the applied R-ratio (Rc omp 
= 0.1). Thus, the value of Rc omp does not represent the actual stress ratio in the matrix and 
the local R m values listed in Table 2 should be used in the analysis. 

Calculation of AK m 

The stress intensity factor range in the matrix, A K m , was calculated using an FB model 
implemented with thermal residual stresses. It is assumed in this model that fiber-matrix 
debonding takes place as a matrix crack progresses past the fibers. The intact fibers in the 
wake of the matrix crack are idealized as a crack bridging pressure, Ap. An unknown but 
constant frictional shear stress, t, is assumed to act along the debond length of the bridging 
fibers. A continuum fracture mechanics analysis and a micromechanical analysis are combined 
using discrete-continuum relations to obtain stress intensity factor solutions for fiber bridging 
matrix cracks of arbitrary size. 

In the calculation of A K m , the crack opening displacement in terms of the bridging 
pressure obtained from the micromechanical analysis and the continuum fracture mechanics 
analysis are combined using discrete-continuum relations to yield a single non-linear integral 
equation in terms of the bridging pressure. Solving for the bridging pressure, AK m can be 
determined. Assuming the discrete-continuum relations in [15], the stress-intensity factor in 
the matrix is: 

A K m = ASjnaF opp - 2 I ApG(x,a,R,W)dx (6) 

Jo. 


was calculated using the VISCOPLY 
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In this equation, AS is the applied far-field stress, F app is the applied stress boundary 
correction factor, a g is the initial unbridged crack length (hole radius or notch length). In 
addition, the term G(x,a,R,W ) is the weight function for the bridging stress intensity factor 
range. Details on the calculation of can be found in [2]. 

Fatigue Life 

A procedure to predict fatigue life is proposed in this study based on the axial stress in 
the first intact fiber calculated using the FB model. The axial stress in the fibers is 
approximated in terms of the bridging pressure (Ap) and the fiber volume fraction (vy): 



( 7 ) 


This equation is valid for fibers closely spaced relative to the crack length. The bridging 
pressure is calculated using the aforementioned fiber bridging model. When the stress in the 
first intact fiber reaches a critical value, failure of that fiber occurs. During fatigue loading of 
TMC, catastrophic fracture of a specimen typically occurs shortly after the onset of fiber 
breakage. By determining the stress in the first intact fiber, the conditions under which a 
specimen will fail, namely, the loading conditions (applied load and number of cycles) and the 
specimen configuration (hole diameter and crack length) can be predicted. 

For this purpose, the onset of fiber breakage was monitored and recorded in several 
cross-ply SCS-6/Timetal-21S specimens containing center holes using the AE technique. A 
detailed description of the AE test set-up and parameters used is provided in [8]. Briefly, AE 
events having high amplitudes (i.e. greater than or equal to 99 dB) were associated with fiber 
breakage. 

The onset of fiber breakage was also determined from the crack growth characteristics. 
A schematic of the crack length as a function of number of cycles is shown in Figure 7. For 
this example, the crack growth is characterized by three regions: (1) initiation and growth, (2) 
slow growth/crack arrest, and; (3) rapid crack growth. After the initial region, the rate of 
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crack growth decreases due to the fiber bridging phenomena. The number of fibers bridging 
the crack increases as the crack length increases thereby reducing the stress-intensity factor in 
the matrix. As the crack extends, the stress in the bridging fibers increases. Under certain 
conditions, the stresses in the fibers may become sufficiently high to cause fiber breakage to 
occur resulting in the rapid crack growth in the third region. In this study, the onset of fiber 
breakage is assumed to occur at the "knee" of the crack growth curve between the second and 
third regions, Figure 7. 

Using the AE technique and the characteristics of the crack growth curve, the 
conditions (i.e.: specimen geometry, crack length, applied load, number of cycles, etc.) under 
which fiber breakage occurs was determined for several specimens. Using these conditions as 
input parameters, the failure stress range for the first intact fiber, Ao^ , was calculated using 

the FB model and was used as a failure criterion to predict fatigue life. These failure stresses 
are the range values ( Ac y 1 — oy* - oy n ) between and S and are listed in Table 3. 

Post -Fatigue Residual Strength 

The post-fatigue residual strength was predicted using two approaches: (1) a fiber 

bridging model, and; (2) a three dimensional finite element program, PAFAC [17]. The axial 
stress in the 0° fiber next to the notch (first intact fiber) was predicted as a function of applied 
load. The composite strength was assumed to be the applied load at which the axial stress in 
the 0° fiber next to the hole reached the fiber strength. 

Selection of the approach used to predict the post-fatigue residual strength depends on 
the characteristics of the damage process zone near the notch. For damage classified as self- 
similar, mode I matrix crack growth, the fiber bridging model is suitable. The FB models are 
applicable to the analysis of SCS-6/Timetal-21S laminates where the damage is primarily a 
single mode I matrix crack. The fiber bridging analysis was previously described. The 
critical fiber strength, af‘, used to determine the onset of fiber breakage in SCS-6/Timetal- 

21S during static loading is listed in Table 3. These failure stresses are absolute values and are 
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lower than the ranges values obtained under fatigue loading. Under fatigue loading, the fiber 
stresses at the minimum applied load, crj"’, is compressive due to the thermal residual stresses. 

Consequently, the fiber stress range under fatigue loading, Ao f “' , is higher than the absolute 
fiber stress under static loading, of 1 . 

In TMC where extensive damage process zones develop, such as in SCS-6/Ti-15-3, a 
more robust approach must be used such as a three-dimensional finite element analysis. 
Details of this approach are outlined in [3-5]. Briefly, a three-dimensional finite element code 
PAFAC [17], is used to calculate the fiber stresses. The damage that evolves during fatigue 
loading, namely, matrix cracking and fiber-matrix debonding, are analytically simulated by 
using reduced stiffness properties of the matrix material. 

RESULTS AND DISCUSSION 

The aforementioned fatigue crack growth methodology was used to determine the 
effects of applied stress, number of fatigue cycles, applied temperature and specimen geometry 
(e.g. stress concentrations and notches) on damage initiation, damage growth, fatigue life and 
residual strength in TMC as discussed in the following sections. 

Matrix Crack Initiation 

Using Equation (1), the effective strain parameter versus the number of cycles to 
matrix crack initiation is plotted in Figure 8 for the SCS-6/Ti-15-3 laminates tested in [1,3]. 
In this figure, CH and DEN designate center hole and double edge notched specimens, 
respectively. The solid line is the baseline curve for the matrix material (same as that shown 
in shown Figure 2). Experimentally, matrix crack initiation was assumed when the crack 
length reached between 0.3-mm and 0.5-mm. As shown in Figure 8, by using Equation (1), 
the experimental data are consolidated along the baseline curve indicating the effective strain 
parameter can be used to predict crack initiation in the matrix reasonably well. 
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Matrix Crack Growth 


The matrix crack growth behavior observed in the two TMC systems was analyzed 
using the FB model implemented with thermal residual stresses [2]. The frictional shear stress 
term, T, in the FB model was used as a curve-fitting parameter to matrix crack growth data, 
reduced using the crack growth law defined by Equation (5). Results are presented in the 
subsequent sections. The effects of thermal residual stresses on the values of r used in fiber 
bridging analysis are also discussed. 

Calculated Crack Growth Curves 

The crack length as a function of number of cycles was calculated and plotted in 
Figures 9 through 12 using values of r which best-fit the experimental results for the two 
material systems having various external load and geometry conditions. The discrete- 
continuum relations developed by Marshall, Cox and Evans (MCE) [13] were used in the 
calculations. As shown for each figure, a good agreement was obtained between the 
experiments (open symbols) and calculations (filled symbols) using a single value of the 
frictional shear stress. 

Of particular interest is the results presented in Figure 12 showing the crack length as a 
function of the number of cycles at room temperature and at an elevated temperature of 200°C 
for two center notched specimens of [0/90] s SCS-6/Timetal-21S. As shown in the figure, 

good agreement was obtained between experiments (open symbols) and the calculations (filled 
symbols) using the values of T = 50 MPa and 15 MPa for the room temperature and elevated 
temperature test results, respectively. The difference in the calculated crack growth behavior 
at the two temperatures is due primarily to the differences in thermal residual stress states 
which effects the calculated values of R m and A K m . 

If the mechanical properties of the matrix are assumed to be identical at both 
temperatures, only R m and A K m affect the crack growth rate defined by Equation (5). At 
elevated temperatures, the axial thermal residual stress in the matrix d r ^ rm is lower than at 
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room temperature which results in lower values of R„. For the room and elevated temperature 
tests, the calculated local R-ratios in the matrix were R m = 0.75 and R m = 0.64, respectively. 
From Equation (5), for a given value of AK m , the crack growth rate increases with an increase 
in the value of R m . Based solely on the values of R m for the two test conditions, the crack 
growth rate should have been higher for the room temperature test than the elevated 
temperature test. Results shown in Figure 12, however, show the opposite trend; the crack 
growth rate data was higher for the elevated temperature test compared to the room 
temperature test. 

Another factor affecting the crack growth rate is the stress-intensity factor range AK m . 
Thermal residual stresses affect the value of A K m . The lower value of t for the elevated 
temperature test translated into a higher value in the stress-intensity factor range as shown in 
Figure 13. At the elevated temperature, the value of AK m was higher than at room 
temperature. Using the crack growth law defined by Equation (1), for a given value of R m , 
the crack growth rates increases with an increase in the value of A K m which agrees with the 
results shown in Figure 12. As a result of the lower thermal residual stresses at the elevated 
temperature test, the value of A K m increased while the value of R„ decreased. The lower 
values of T for the elevated temperature results suggests that the radial thermal residual 
stresses acting across the debonded fiber-matrix interfaces in the 0° plies were lower. These 
radial stresses "choke" the fibers which enhances the shear transfer between the fiber and 
matrix. By lowering the radial thermal stresses, the amount of shear transfer and displacement 
constraint between the fiber and matrix were lowered thereby increasing the stress-intensity 
factor. The increase in the value of AK m had a more profound effect on the crack growth rate 
than the decrease in the value of R„. 

Effect of Thermal Residual Stress on the Frictional Shear Stress 

The effect of thermal residual stresses on fiber bridging analysis is revealed in the 
values of the frictional shear stress term required to reduce the crack growth data which are 
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listed in Table 2. The values of T listed in this table were determined using the discrete- 
continuum relations developed by McMeeking and Evans (ME) [15] with and without thermal 
residual stress. For the unidirectional SCS-6/Ti-15-3 specimens tested in [6], the value of T 
required to fit crack growth data varied from 12 MPa to 41 MPa when thermal residual 
stresses are neglect (i.e. original fiber bridging analysis). On the other hand, when thermal 
residual stresses were included in the analysis, only a single value of frictional shear stress, T 
= 60 MPa, was required to curve-fit the data. Thus, the scatter in the value of r was reduced 
significantly by considering thermal residual stresses in the fiber bridging analysis. In fact, for 
a given material system, lay-up and temperature, a single value of T was sufficient to analyze 
the crack growth data. Thermal residual stresses are an important factor overlooked in the 
original formulations of the FB models. 

Fatigue Life 

Using the values of frictional shear stress used to evaluate the crack growth data 
presented in the pervious section, the fatigue life was estimated based on the stress in the first 
intact fiber. The critical value of stress in the fiber at which fiber breakage occurred was first 
determined for several specimens of cross-ply SCS-6/Timetal-21S. For these specimens, the 
Acoustic Emission (AE) technique and the characteristics of the crack growth curve were used 
to determine the conditions (i.e.: specimen geometry, crack length, applied load, number of 
cycles, etc.) under which fiber breakage occurred as described previously. Using these 
conditions as input parameters, the value of stress in the first intact fiber when fiber breakage 
occurs was calculated for both static and fatigue applied loadings using the FB model. Results 
are presented in Table 3. As shown, the value of the axial stress in the first intact fiber was 
nearly a constant value for the specimens listed. In this study, the average value of the stresses 
listed in Table 3 was used as a critical value to determine the onset of fiber breakage. Under 
fatigue loading conditions, the value of stress used was 2500 MPa and under static loading 
conditions, the value used was 2000 MPa. 
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Having the critical value of fiber stress, the life prediction methodology was applied to 
a [0/90 ] 2s SCS-6/Timetal-21S laminate containing a center hole subjected to S max = 200 MPa 
and Rcomp = 0.1. The crack length as a function of the number of cycles was monitored and 
is plotted in Figure 1 1 along with the calculation made using the FB model and a value of t = 
50 MPa. Crack growth was recorded up to 200,000 cycles. This specimen failed after 
222,600 cycles. The stress in the first intact fiber as a function of the number of fatigue cycles 
was calculated and plotted in Figure 14. As shown in the figure, the stress in the first intact 
fiber reaches the average critical value (2500 MPa) after 285,000 cycles. 

Post -Fatigue Residual Strength 

The post-fatigue residual strength was evaluated in two cross-ply TMC laminates 

containing center holes, namely, in a SCS-Ti-15-3 specimen and a SCS-6/Timetal-21S 
specimen. The SCS-6/Ti-15-3 specimen had been subjected to S mgx = 250 MPa for 200,000 

cycles during which extensive matrix cracking occurred [3]. This specimen was loaded quasi- 

statically in tension to failure and a residual strength of 325 MPa was measured. 

The SCS-6/Timetal-21S specimen had been subjected to S ma = 150 MPa for 200,000 

cycles during which two symmetric cracks grew from the hole. This specimen was then 
loaded quasi-statically in tension to failure and a residual strength of 350 MPa was measured. 
The AE technique was used to determine the occurrence of fiber breakage. For this specimen, 
the number of high amplitude events, that is, AE events having an amplitude greater than 
98dB, associated with fiber breakage, is plotted in Figure 15 as a function of the applied load. 
As shown, fiber breakage occurs at an applied load of 180 MPa. The number of high 
amplitude events increases up to failure at 350 MPa. This indicates that stable fiber breakage 
occurred prior to catastrophic fracture. From the onset of fiber breakage, an increase in 
applied load was required for subsequent fiber breakage's. 
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Fiber Bridging Analysis 

The FB model was used to evaluate the post-fatigue residual strength for SCS- 
6/Timetal-21S specimen. Since the fatigue damage that occurred in SCS-6/Timetal-21S 
specimen can be considered as mode I self-similar crack growth, the application of the FB 
model in the analysis was suitable. The axial stress in the first intact fiber as a function of the 
applied load was calculated. The onset of set of fiber breakage was assumed to occur when 
the stress in the first intact fiber equaled the critical value of 2000 MPa. The critical value 
was calculated as described previously. Using the FB model, the stress in the first intact fiber 
reaches the critical value at an applied load of 235 MPa. Using the AE technique, the 
emission of high amplitude events associated with fiber breakage occurred at 180 MPa. A 
comparison is shown in Figure 15. While both these values are in agreement, this specimen 
failed catastrophically at 350 MPa. This indicates that the occurrence of first fiber breakage 
may not be a good indicator of post-fatigue residual strength for the SCS-6/Timetal-21S. The 
results shown in Figure 15 suggest that a stable fiber breakage occurred prior to catastrophic 
fracture. 

PAFAC Analysis 

A three dimensional finite element analysis (PAFAC) was used to predict the virgin 
static and post-fatigue residual strengths in cross-ply SCS-6/Ti-15-3 laminates containing 
center holes. For the virgin specimen, the constituent properties shown in Table 1 were used. 
In addition, fiber-matrix debonding in the 90° plies was modeled as described in [5]. The 
damage that occurred in fatigued specimen was quite extensive and widespread consisting of 
multiple matrix cracks and fiber-matrix debonding. This damage was simulated using a 
reduced matrix modulus as described in [3]. 

The axial stress in the 0° fiber next to the center hole as a function of the applied stress 
is plotted in Figure 16 for both cases. The horizontal line indicates the assumed fiber strength 
of 4400 MPa, and the two vertical lines show the experimental strengths of the two conditions. 
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The fiber strength was calculated from the measured strain to failure of an [0/90] s unnotched 
specimen (ej!^ = 0.0011) multiplied by the fiber modulus (400 GPa). As shown in Figure 

16, the applied stress at which the 0° fiber stress equals the fiber strength corresponds closely 
to the measured strength for both conditions. Using the axial stress in the 0° fiber next to the 
center hole as the criterion for laminate failure, the predicted strengths for the virgin and post- 
fatigue conditions were 483 MPa and 301 MPa, respectively. Both predictions are within 8% 
of the experimental values. 

SUMMARY 

In this study, a micromechanics-based life prediction methodology was presented and 
used to determine the effects of applied stress, number of fatigue cycles and specimen 
geometry on damage initiation, damage growth and residual strength in notched Titanium 
Matrix Composite (TMC). This methodology was applied to two TMC systems, namely, 
SCS-6/Ti-15-3 and SCS-6/Timetal-21S. 

Damage initiation was predicted using a local effective strain approach. A finite 
element analysis of a SCS-6/Timetal-21S monolayer containing a center hole verified the 
prevailing assumptions made in the formulation of this model. The finite element analysis 
showed that the primary stress concentration was due to the center hole and that the local stress 
concentrations due to machined fiber ends were secondary. Using the effective strain 
approach, the conditions under which matrix cracks initiated were predicted. 

Damage growth analysis of fiber bridged matrix cracks was developed based on a fiber 
bridging (FB) model which accounts for thermal residual stresses. The FB model combines 
continuum fracture mechanics and micromechanics analyses to derive stress intensity factor 
solutions of matrix cracks bridged by fibers. The unknown constant frictional shear stress 
assumed in this model was used in this study as a curve fitting parameter to the available 
experimental data. For a given material system, lay-up and applied temperature, a single 
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value of frictional shear stress was sufficient to accurately model the matrix crack growth 
behavior. 

The fatigue life was predicted based on the axial stress in the first intact 0° fiber 
calculated using the FB model. The critical value of stress at which fiber breakage occurred 
was first determined for several specimens. For these specimens, the Acoustic Emission (AE) 
technique and the characteristics of the crack growth curve were used to determine the 
conditions under which fiber breakage occurred. Using these conditions as input parameters, 
the value of stress in the first intact fiber when fiber breakage occurs was calculated using the 
FB model and was found to be nearly a constant value. This critical value was then used to 
predict the onset of fiber breakage during fatigue loading in a subsequent specimen. It was 
assumed that fatigue life is given by the number of fatigue cycles at which fiber breakage 
initiates. There was approximately a 22% difference between the observed and predicted 
fatigue lives. 

The post-fatigue residual strength was predicted based on the axial stress in the first 
intact 0° fiber calculated using either the FB model and a three-dimensional finite element 
analysis, depending on the classification of fatigue damage. The fatigue damage developed in 
the SCS-6/Timetal-2 1 S laminates can be idealized as a mode I, self-similar crack. For this 
case, the fracture mechanics-based approaches are appropriate and the FB model was used. 
Using the FB model, the predicted applied load causing the onset of fiber breakage was in 
agreement with AE results. There was approximately a 23% difference in the results. 

The virgin static and post-fatigue residual strength was predicted for cross-ply SCS- 
6/Ti-15-3 laminates containing a center hole using a three-dimensional finite element analysis. 
The fatigue damage that occurred was quite extensive and widespread consisting of multiple 
matrix cracks and fiber-matrix debonding. This damage was simulated using a reduced matrix 
modulus. Using the axial stress in the 0° fiber next to the center hole as the criterion for 
laminate failure, the predicted strengths were within 8% of the experimental values. 
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Table 1. Material Properties: 


Material System 


Property 

Matrix E m (GPa) 

v* 

C 

n 

m 

a m mm/mm/C 

Fiber: E f (GPa) 

v / 

a f mm/mm/C 

Composite: Lay-Up 

v / 

E l (GPa) a 
E T (GPa) a 
G lt (GPa) a 
v a 

G r * trnl (MPa) a 

Lay-Up 

v / 

E l (GPa) a 
E t (GPa) a 
G lt (GPa) a 
v a 

V LT 

d^ ,rm (MPa) a 


SCS-6/Ti-15-3 

SCS-6/Timetal-21S 

92.4 

116.0 

0.34 

0.34 

3.8E-11 

5.5E-11 

2.8 

2.5 

0.35 

0.35 

9.7x10-* 

8.6 x 10" 6 

400.0 

400.0 

0.25 

0.25 

3.6 x 10"* 

3.6x10-* 

[0] 8 

[0/90] s 

0.33 

0.36 

192.2 

188.4 

130.1 

188.4 

45.9 

58.1 

0.32 

0.26 

202.0 b 

221.9 c 

[0/90] s 

[0/90] 2s 

0.35 

0.38 

168.4 

192.9 

168.4 

192.9 

47.4 

59.3 

0.25 

0.26 

183.8 b 

230.5 c 


Predicted using VISCOPLY [9] based on constituent properties 

Axial .thermal residual stress in matrix in 0° plies for cool-down from 555 °C 

Axial .thermal residual stress in matrix in 0° plies for cool-down from 620 °C 
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Table 2. Specimen loading history and fiber bridging analysis parameters: 
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Frictional shear stress required to analyze data using ME relations [15] when thermal residual stresses are neglected 
Not applicable 

Hole diameter to width ratio (d/W) 

Frictional shear stress required to analyze data using ME relations [15] neglecting the effect of parallel cracks 


Table 3. Calculated Stress in First Intact Fiber at the Onset of Fiber Breakage 
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Static Strength Using 
Finite Element Method 
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Figure 1. Micromechanics-based approaches used to predict and evaluate damage initiation, damage growth, 
and residual strength in Titanium Matrix Composites (TMC). 
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Figure 2. Effective strain parameter as a function of the number of cycles to failure for as-received Ti-15-3 matrix material [1], 
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Figure 3. Schematics of cracking pattern 
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ises stress distribution in vicinity of machined fiber ends. Maximum stress and strain concentration is 
center hole. Machined fiber ends cause secondary stress and strain concentrations 




CONTOUR LEVEL 
(MPa) 


in 

o> 

CO 






ID 

CM 

N 


in 

o 

T- 

O 

o 

K 

CO 

CO 

CO 

T* 

T " 


O 

o> 

CO 

00 

II 

II 

II 

1! 

II 

II 

II 

< 

CD 

O 

O 

LU 

Urn 

O 


CO 

N 

«|— 

CD 

o 

** 

00 

CM 

o> 


o 

in 

T" 

CD 


h» 

N 

N 


CD 

CD 

in 

in 

rr 

II 

II 

II 

II 

II 

II 

ii 

ii 

x 

mmm 

"D 


-1 

s 

Z 

O 



ip behind the second cut 




Failure 


fo 

o U) 


.*= (0 
O ^ 



33 


Figure 7. Schematics of typical crack growth curve observed in TMC. Characterized by three regions: 1). Initiation 
and rapid crack growth; 2. Slow crack growth, and ; 3. Rapid crack growth up to failure. 
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Number of Cycles 

Figure 9. Crack growth data for center notched [0] g SCS-6/Ti-15-3 laminates at several applied stress levels. 
Comparison between experiments and calculations made using Fiber Bridging Model. 
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Figure 10. Crack growth data for center hole [0/90] s SCS-6/Ti-15-3 laminates at several applied stress levels. 
Comparison between experiments and calculations made using Fiber Bridging Model. 





37 


Figure 11. Crack growth data for center hole cross-ply SCS-6/Timetal-21S laminates at several applied stress levels. 
Comparison between experiments and calculations made using Fiber Bridging Model. 
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Figure 12. Crack growth data for center notched [0/9Q] SCS-6/Timetal-21S laminates at two temperature levels. 
Comparison between experiments and calculations made using Fiber Bridging Model. 
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laminates at two temperature levels. 
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Figure 14. Fatigue life prediction based on stress in the first intact fiber calculated using the fiber 
bridging model. 
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Figure 15. Determination of fiber breakage via Acoustic Emission technique and fiber bridging model. 
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Figure 16. Axial stress in 0° fiber next to hole as a function of the applied stress for virgin and post-fatigue conditions in 
[0/90] s SCS-6/Ti-15-3 laminate, vf = 0.355. Applied stress at which 0° fiber stress reaches the fiber strength 
corresponds to the notched strength for both conditions. 
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